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Abstract. High-energy neutrinos could be produced in the interaction of charged cosmic rays with 
matter or radiation surrounding astrophysical sources. Even with the recent detection of extrater¬ 
restrial high-energy neutrinos by the IceCube experiment, no astrophysical neutrino source has yet 
been discovered. Transient sources, such as gamma-ray bursts, core-collapse supernovae, or active 
galactic nuclei are promising candidates. Multi-messenger programs offer a unique opportunity to 
detect these transient sources. By combining the information provided by the ANTARES neutrino 
telescope with information coming from other observatories, the probability of detecting a source is 
enhanced, allowing the possibility of identifying a neutrino progenitor from a single detected event. 
A method based on optical and X-ray follow-ups of high-energy neutrino alerts has been developed 
within the ANTARES collaboration. This program, denoted as TAToO, triggers a network of robotic 
optical telescopes (TAROT and ROTSE) and the 5'w//i-XRT with a delay of only a few seconds after 
a neutrino detection, and is therefore well-suited to search for fast transient sources. To identify an 
optical or X-ray counterpart to a neutrino signal, the images provided by the follow-up observations 
are analysed with dedicated pipelines. A total of 42 alerts with optical and 7 alerts with X-ray im¬ 
ages taken with a maximum delay of 24 hours after the neutrino trigger have been analysed. No 


optical or X-ray counterparts associated to the neutrino triggers have been found, and upper limits on 
transient source magnitudes have been derived. The probability to reject the gamma-ray burst origin 
hypothesis has been computed for each alert. 

Keywords: ANTARES - Neutrino astronomy - Transient sources - Gamma-ray bursts - Optical/X- 
ray follow-up 

1 Introduction 

High-energy neutrinos are expected to be produced in the interactions of accelerated charged cosmic 
rays with matter and radiation fields within and surrounding their astrophysical sources [1,2]. Neu¬ 
trinos are unique messengers for studying the high-energy Universe as they are neutral and stable, 
interact weakly, and travel directly from their source without absorption. Neutrinos easily escape the 
acceleration regions and propagate through space without interaction, unlike charged particles, which 
are deflected by magnetic fields and no longer point back to their source. Detection of high-energy 
neutrinos from an astrophysical source would be a direct signature for the presence of hadronic ac¬ 
celeration, and would therefore provide crucial information on the origin of very high-energy cosmic 
rays. 

The production of high-energy neutrinos has been proposed to occur in several kinds of as¬ 
trophysical sources, in which the acceleration of hadrons may occur, such as active galactic nuclei 
(AGN), gamma-ray bursts (GRBs), supernova remnants and microquasars. Many of these astro- 
physical accelerators show transient behaviour. Variations in the energy output of the most powerful 
astrophysical objects cover a large range in the time domain, from seconds for GRBs [3] to weeks for 
AGN [4] or core-collapse supemovae (CCSNe) [5]. The particularity of these high-energy phenom¬ 
ena is that they radiate over the entire electromagnetic (EM) spectrum, from the radio domain to TeV 
gamma rays. Other approaches to astronomy do not use EM radiation at all however, but use other 
messengers, such as neutrinos, gravitational waves or cosmic rays. The detection of astrophysical 
sources of these other messengers is very difficult due to the very small number of expected events 
and a large background contamination. A way to overcome this difficulty is to combine the detection 
of non-EM messengers with the EM signal, to provide a multi-messenger dataset. 

The two largest operating large-volume neutrino telescopes in the world are: the Northern hemi¬ 
sphere ANTARES telescope [6], located in the Mediterranean Sea, and the Southern hemisphere 
IceCube telescope [7], located in the South Pole ice. These telescopes are designed to search for 
high-energy cosmic neutrinos (Ey > 100 GeV) generated in extreme astrophysical sources. Each 
detector is able to monitor a full hemisphere of the sky (or even the whole sky if downgoing events 
are considered) with a high duty cycle. The recent IceCube discovery of extraterrestrial high-energy 
neutrinos [8, 9] opened new windows in the field of astroparticle physics. This result has fixed the 
scale of neutrino fluxes in the Universe and has inspired a large number of hypotheses for their origin, 
mainly due to the poor localisation of the events. Up to now, no high-energy neutrino source has been 
identified. 

Searches for transient astrophysical phenomena offer very promising opportunities for high- 
energy neutrino telescopes, because the relatively short duration of the events means that the level of 
contamination (from atmospheric muons and atmospheric neutrinos) is strongly reduced. Targeted 
searches performed by the ANTARES telescope on GRBs [10], on AGN flares [11] and microquasar 
outbursts [12] have so-far yielded limits on the neutrino production in these sources. Taking full 
advantage of the possibilities offered by multi-messenger searches for transient sources, a multi¬ 
wavelength follow-up program, denoted as TAToO (Telescopes-ANTARES Target of Opportunity), 
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has operated within the ANTARES Collaboration since 2009 [13]. It is based on optical and, since 
mid 2013, X-ray follow-ups of selected high-energy neutrino events very shortly after their detection. 
The network is composed of the small robotic optical telescopes TAROT [14] and ROTSE [15], and 
the Swift X-ray telescope (XRT) [16]. This approach has the advantage that it does not require an a 
priori hypothesis on the nature of the underlying neutrino source. To be sensitive to all types of time 
variability in the astrophysical sources, the observational strategy is composed of a real-time obser¬ 
vation for rapidly fading sources, such as GRB afterglows, complemented by several observations 
during the following month specially adapted to detect the rising light curve of CCSNe. Unprece¬ 
dented in this domain, the ANTARES telescope is able to generate alerts within a few seconds after 
the neutrino detection, and provides a precision on the reconstructed direction to better than 0.5° at 
high energies (E > 1 TeV). The TAToO system is therefore well-suited to searching for rapid tran¬ 
sients showing time variability at the minute scale. A similar program is running in IceCube [17] 
since 2008, but with larger alert generatoin delay. 

With this type of real-time analysis, one neutrino associated with a transient optical or X-ray 
counterpart would be a significant event. Since 2009, around 150 neutrino alerts have been sent to 
the optical telescope network. Among them, 42 alerts have had follow-up observations within one 
day after the neutrino detection. With the Swift-XRT follow-up, a total of 7 alerts have been studied. 

In this paper, the first results of the analysis of the rapid follow-up observations associated with 
the TAToO neutrino alerts are presented. A future dedicated paper will discuss the results of the 
analysis of the long-term follow-up. A brief description of the ANTARES experiment and the alert 
system is given in section 2. The neutrino data set used in this analysis is described in section 3. 
The optical and X-ray analyses, as well as the results are presented in sections 4 and 5 respectively. 
Einally, the implications of these results for GRB neutrino production models are discussed in section 
6 . 

2 TAToO and its multimessenger network 
2.1 The ANTARES neutrino telescope 

Eor the ANTARES telescope [6], events are detected underwater by the Cherenkov light induced in 
the medium by secondary particles resulting from neutrino interactions in the darkness of the deep 
sea. The detector is a 3D array of roughly 900 optical sensors in an approximately cylindrical volume 
of base diameter 190 m and height 350 m. Each optical sensor includes a 10-inch photomultiplier 
oriented downwards to maximise the detection efficiency to upgoing particles. 

Owing to their low interaction probability, neutrinos have the ability to cross the Earth, and 
muon neutrinos interacting close to the detector produce an upgoing muon. At high energy (Ey > 1 
TeV), the Eorentz boost is such that the muon and neutrino directions are identical to within a few 
tenths of a degree, so that accurately reconstructing the muon trajectory provides a good estimate 
of the arrival direction of the neutrino candidate. Another source of muons is due to cosmic rays 
that hit atmospheric nuclei and produce particle showers. These ’’atmospheric muons” have down¬ 
going trajectories, which can therefore be used to discriminate against them. These muons, whose 
abundance at the ANTARES detector is roughly six orders of magnitude larger than those induced 
by atmospheric neutrinos, are the main background and have to be efficiently suppressed. The signal 
recorded by the detector is also polluted by decays of natural Potassium 40 in sea water and biolumi¬ 
nescence [6]. These backgrounds can be mitigated by applying causality criteria to the observed light 
patterns. 

The telescope, fully equipped since May 2008, is anchored to the seabed at 2475 m below the 
Mediterranean sea surface, at 42°48’N 6°10’E, with a maximum efficiency in the Southern sky [18]. 
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All signal hits from the photomultipliers are digitised in situ with a time resolution of about 1 ns, 
time stamped by a clock system synchronous to a standard GPS and sent to a 40 km distant shore 
station. When several optical sensors show a topology compatible with a muon track crossing the 
instrumented volume, all photomultiplier data is stored by the PC farm at the shore station and made 
available within a few seconds for further track reconstruction and physics analysis. 

2.2 The TAToO alert system 

The criteria for the TAToO trigger are based on the features expected from astrophysical sources. Sev¬ 
eral models predict the production of neutrinos with energy greater than 1 TeV from GRBs [19-22], 
CCSNe [23] and AGN [24] . A basic requirement for the near coincident observation of a neutrino 
and an optical counterpart is that the pointing accuracy of the neutrino telescope should be at least 
comparable to the field of view of the TAROT and ROTSE telescopes (w 2° x 2°). After the selection 
of upgoing events, which removes the huge background of atmospheric muons, the ANTARES neu¬ 
trino sample consists mainly of atmospheric neutrinos. Different criteria are used to select candidates 
with an increased probability to be of cosmic origin [13]. A fast and robust algorithm is used to 
reconstruct the data [25], which uses an idealised detector geometry and is thus independent of the 
dynamical positioning calibration. This reconstruction, and a subsequent quality selection (Nline> 2, 
zenith angle < 0, and Q < 1.3-1- [0.04(Ahits - 5)]^, see Ageron et al. 2012) allows the rate of events 
to be reduced from few Hz down to few mHz. The remaining events are then passed to a more pre¬ 
cise reconstruction tool [26] which allows the neutrino nature of the event to be confirmed, and the 
angular resolution to be improved. 

Three online neutrino trigger criteria are currently implemented in the TAToO alert system: 

• Doublet trigger: the detection of at least two neutrino-induced muons coming from similar 
directions (< 3°) within a predefined time window (< 15 minutes). 

• High-energy trigger: the detection of a single high-energy (~ 7 TeV) neutrino-induced muon. 

• Directional trigger: the detection of a single neutrino-induced muon for which the direction 
points toward a local galaxy (< 0.5°). The coordinates of the galaxies are selected from the 
GWGC catalogue [27] with a distance cut at 20 Mpc. 

In agreement with the optical telescopes, the total trigger rate is tuned to 25 per year. This 
rate is dominated by high-energy and directional triggers, as until now no doublet trigger has been 
sent to the network. The accidental coincidence rate due to two uncorrelated events is estimated to 
be 7 X 10~^ per year. The high-energy trigger typically requires more than 70 photomulitplier hits 
and a total amplitude greater than 150 photoelectrons. To comply with the rate of 6 alerts per year 
required by the Swift satellite, a subset of the high-energy trigger, denoted as the very high-energy 
trigger, provides a dedicated trigger for the XRT. It typically requires more than 80 hits and 300 
photoelectrons. The bidimensional distribution of variables used for the selection of both high-energy 
and very high-energy triggers is illustrated in figure 1 for neutrino candidates recorded from 2012 to 
March 2015. Neutrinos falling outside the red box can lead to an alert if they fulfil the directional or 
the doublet trigger. The directional trigger was implemented in late 2011. The performances of these 
three triggers are described in table 1 . 

Eigure 2 shows the estimate of the point spread function (PSE) for a typical high-energy neutrino 
alert, compared to the field of view of TAROT and ROTSE. Eor the highest energy events, the angular 
resolution reaches less than 0.3° (median value). 

Before 2012, the alert system sent alerts about 50 seconds after the neutrino events. After 2012, 
a major DAQ system improvement enabled ANTARES to send alerts a few seconds (~ 3 - 5 seconds) 
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Number of PMT hits 


Figure 1. Scatter plot of the number of photomultiplier hits vs. the total amplitude (in photoelectrons) for 
neutrino candidates. The high-energy trigger criterion is fulfilled by events inside the red box. Events inside 
the blue box correspond to the very high-energy trigger for the Swift-XRT. 


Table 1. Performances of the three alert criteria. 


Trigger 

Angular Resolution (median) 

PSF coverage” 

Atmospheric muon contamination 

Mean energy'’ 

High energy 

0.25 - 0.3° 

96% 

< 0.1 % 

~7TeV 

Directional 

0 

O 

1 

o 

90% 

~ 2 % 

~ 1 TeV 

Doublet 

< 0.7° 


0% 

~ 100 GeV 


“ Percentage of the ANTARES PSF covered by a 2° x 2° telescope field of view, assuming events produced in GRBs. 
Neutrino energy weighted assuming the atmospheric muon neutrino spectrum. 


after the detection of the neutrinos. Figure 3 displays this latency for the 150 alerts collected since 
the commisioning of the TAToO alert system. 


3 Neutrino data set 

This dataset is composed of 42 alerts with an early optical follow-up and 7 alerts with an X-ray follow¬ 
up. Only one alert has an early follow-up with both optical and X-ray telescopes. The direction, time 
and energy of the alerts are reported in table 2. A total of 34 and 14 neutrinos have been triggered 
as high-energy and directional alerts respectively. These numbers are in agreement with expectations 
from Monte Carlo simulation assuming that they are of atmospheric origin. 

A skymap in galactic coordinates of the neutrino alert directions which have triggered the alert 
system since 2009 is displayed in figure 4. 
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(deg) 

Figure 2. Bi-dimensional angular resolution for a typical high-energy neutrino alert. The black square cori'e- 
sponds to the TAROT and ROTSE telescope field of view (=s 2° x 2°). 



Figure 3. Latency for the 150 triggers sent by the alert system. The step from 50 to 5 seconds corresponds to 
an upgrade of the ANTARES DAQ system. 


4 Optical follow-up 

4.1 The optical telescope network 

The optical follow-up of ANTARES alerts is performed using ground-based telescopes dedicated to 
early observations of GRBs. Six such telescopes are involved: ROTSE-III (USA, Namibia, Australia, 
Turkey) and TAROT (Erance and Chile). ROTSE-III is a network of four identical 0.45 m tele- 
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Figure 4. Skymap in galactic coordinates of the candidates which triggered the alert system. Black triangles 
refer to the triggers which had an early optical follow-up, while red triangles show the triggers with no early 
observation. The visibility of the ANTARES detector is represented in greyscale, with a visibility of 100 % for 
the darkest zone and a white region with no visibility. 


scopes [15] and TAROT is a network of two identical 0.25 m telescopes [14]. ROTSE has stopped 
its activity progressively over the last years: ROTSE 3a (Australia), 3b (Texas), 3c (Namibia) and 
3d (Turkey) have been stopped in 07/2011, 09/2014, 12/2010 and 12/2012, respectively. These tele¬ 
scopes had a field of view of 1.9° x 1.9° and a spatial sampling of 3.3 arcsec/pixel. The sensitivity 
of the TAROT and ROTSE telescopes is about the same: for a signal to noise ratio of 5, an exposure 
time of 180 seconds and a clear filter, the limiting magnitude is ~ 18.5. 

Alerts are sent automatically by an ANTARES server using a socket protocol, similar to the 
Gamma-ray bursts Coordinates Network (GCN) system developed by [28] for GRBs. If we consider 
that ANTARES events originate from GRBs, the rapid decay of their afterglow light implies the use 
of telescopes with very fast slewing. The time reaction of these optical telescopes is of the order of a 
few seconds when the field of view is observable. Eigure 5 shows a map of fhe real-time visibilify of 
ANTARES neufrino alerts by a single felescope as a funcfion of ifs location on Earfh. 

Eor fhe prompf observafions, fhe efficiency of an opfical ground-based felescope depends on fhe 
objecf visibilify (nighf and field elevafion) and fhe local condifions (weafher and technical problems). 
Eor insfance, among 744 GRBs friggered by Swift befween Sepfember 2006 and January 2015, only 
18 % (i.e. 134) were accessible during fhe prompf emission af TAROT Chile. Only 45 GRBs were 
really observed during fhis period. This implies an efficiency of abouf 33 % for TAROT Chile sife 
(i.e. weafher and fechnical problems). As a comparison, for fhe ANTARES alerfs, figure 5 shows 
fhaf only 20 % of friggers are prompfly accessible from TAROT Chile (i.e. 30 among 150). Table 3 
shows fhaf only 11 ANTARES friggers were really observed in less fhan a few minufes by TAROT 
Chile, indicafing a similar efficiency for ANTARES and Swift friggers. The combinafion of visibilify, 
weafher, and reliabilify of TAROT Chile allows fhe prompf observation of one frigger every 15. This 
emphasizes fhe necessity of using a telescope nefwork fo cover a prompf follow-up of fhe ANTARES 
friggers. 

The early observation sfrafegy of TAROT consisfs of faking 6 images of 180 seconds exposure 
as fasl as possible while for ROTSE, up fo 30 images of 20 fo 60 seconds are faken (depending on 
fhe presence of fhe Moon in fhe observed sky). Each image is preprocessed direcfly on sife and fhen 
fransferred fo fhe ANTARES dafa storage for fhe image analysis. 
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Figure 5. The contours of the world map indicate the percentage of neutrino triggers visible immediately 
(based on 140 ANTARES alerts). Blue and magenta points are the locations of TAROT and ROTSE telescopes 
respectively. The black cross indicates the antipodal point of the ANTARES experiment. 


4.2 Optical image analysis 

To search for a transient optical counterpart to the neutrino trigger, images provided by the follow-up 
of neutrino alerts must be processed and analysed (TAROT and ROTSE images are always analysed 
separately). For this, a new pipeline based on an image subtraction method has been developed 
(figure 6). The first step of this program is the astrometric and photometric calibration of images. 

A necessary step of the image subtraction method is to calibrate and align all the images. To 
calibrate astrometry, SExtractor [29] is used to extract sources from an image and to build a cata¬ 
logue containing the coordinates of each detected source. The extraction threshold is set at 2.5 cr 
above background noise. This catalogue and the information provided by the image header (frame 
dimensions. World Coordinate System data, etc.) are then read by SCAMP [30] and compared with a 
reference catalogue to align the images. 

The photometric calibration is done with LePHARE [31]. This program can compute theoretical 
magnitudes for different bands by using a reference photometric catalogue and a list of filter response 
curves. The photometric catalogue is provided by NOMAD [32] and contains magnitudes in B, V, 
R, J, H, K bands for the brightest stars in the field of the image. Knowing the charge-coupled device 
(CCD) response curve, magnitudes in the CCD band^ can be computed. 

After calibration, one can select a reference image (REF) to perform the subtraction. This 
reference must be a very good quality image (in terms of limiting magnitude) with no or less signal 
than in the image to analyse. As we are looking for transient sources with a high variability on 
the minute/hour timescale, only images (IM) taken less than 24 hours after the neutrino trigger are 
analysed. Therefore, a reference can be chosen among images taken several days or weeks after 
the alert. Some corrections are then applied to images (e.g. gaussian filter, flux normalisation) and 
subtraction IM-REF is made pixel by pixel using the imarith function from the IRAE package [33]. 
Counterpart candidates are automatically selected among positive residual objects on the subtracted 
image if they satisfy one of the two following criteria: 

• New source: the candidate is detected in IM and not in REF. 

• Magnitude variation: the source in IM is at least 0.5 magnitude brighter than in REF. 

'All our images are taken with the clear filter. 


- 7 - 













All these candidates are then checked by eye and light curves are built for the interesting ones. 
Only candidates with at least two points in the light curve (i.e. candidates detected in at least two 
images), or one very bright point (i.e. ~ 2 mag below the limiting magnitude) are considered. Aster¬ 
oids and variable stars are rejected. If a candidate showing variability (for example, a fast decrease 
as typical for GRB afterglows) is not already catalogued, it is considered as a strong candidate to be 
the counterpart. 


Image from follow-up Reference Residue 



t_t 


Alignment, flux normalisation, 
gaussian filtering... 

Figure 6. Optical image analysis based on the subtraction method. A good reference image (center) is sub¬ 
tracted from an image (left) where the counterpart signal is expected. Positive residues in the subtracted image 
(right) are then analysed. 

The detection efficiency of the image analysis pipeline has been estimated by injecting fake 
GRB afterglows in three different observations consisting of 6 images of 180 seconds exposure each. 
The mean limiting magnitude of these images is 16.5, 17.9 and 18.7 for the three observations, re¬ 
spectively. Assuming a light curve decay oc r\ 125 GRB afterglows between magnitude 12 and 19 
have been inserted in each of the three sets of images. The efficiency, presented in figure 7, is given 
by the fraction of inserted GRB afterglows that has been detected by the processing and the candidate 
identification. The efficiency of the pipeline is better than 95 % for bright GRB afterglows (magni¬ 
tude between 12 and 14) and reaches ~ 50 % around 1 magnitude below the limiting magnitude of 
images (vertical dashed lines). 

4.3 Optical follow-up results 

A total of 42 alerts from January 2010 to January 2015 with early images have been analysed. The 
telescopes that observed the regions around the neutrino direction, the number of early images anal¬ 
ysed, their exposure as well as the delay and the limiting magnitude of the first image, are given in 
Table 3. For 11 alerts, the time between the neutrino detection and the start of the acquisition of the 
first image is less than 1.2 minutes and is as low as 25 seconds for 7 alerts. 

No optical counterpart associated with one of the 42 neutrinos was found. This could indicate 
that the observations were made too late when the sources could have faded below the sensitivity 
of the optical telescopes or that the sources occurred outside the field of view of the telescopes, or 
that the neutrinos that generated the triggers have another origin. Upper limits on the magnitude of 
possible transient sources which could have emitted the neutrino have been derived. These limits 
correspond to the limiting magnitude of images that is the faintest signal that can be detected. As we 
are looking for rapidly-fading sources, the signal is supposed to be stronger in the first image of the 
observation, so the upper limits are the limiting magnitude of each first image computed at the 5 cr 
level and corrected for galactic extinction [34] (table 3). 






Figure 7. Detection efficiency of the optical image analysis pipeline estimated by inserting fake GRB after¬ 
glows into three sets of images with limiting magnitudes (vertical dashed lines) of 16.5 (blue), 17.9 (red) and 
18.7 (black) and considering a decay of the light curve oc . 

5 X-ray follow-up 

5.1 The Swift X-ray telescope 

The Swift satellite [35] with its XRT [16] provides a unique opportunity to observe X-ray counterparts 
to neutrino triggers. The detection sensitivity of the XRT is 5 x 10“^^ erg cm“^ s“' in 1 ks, with an 
energy band covering from 0.3 to 10 keV. The follow-up of ANTARES alerts with the Swift-XRT (see 
section 2.2) was implemented in June 2013. Due to the small field of view (radius ~ 0.2°) of the XRT 
and the typical error radius of an ANTARES alert (~ 0.3 - 0.4°), each observation is composed of 4 
tiles up to 2 ks exposure each. This mapping covers about 72 % of the ANTARES PSE for a very high- 
energy neutrino (figure 8). The choice of these 4 tiles corresponds to a trade-off between the requested 
exposure time and the neutrino PSE coverage. With an exposure of 2 ks, a sensitivity of 2.5 x 10“^^ 
erg cm“^ s“^ can be reached [36, figure 14]. Eor each trigger, a single Automatic Target is uploaded, 
and XRT observations start as soon as possible. Swift automatically divides the visibility window of 
the ANTARES trigger location in each spacecraft orbit between the 4 tiles. This is repeated on each 
orbit until the requested exposure time has been gathered for each tile. An automatic analysis of the 
data, described in the following section, is performed soon after, and in the case where an interesting 
candidate counterpart is found, further observations are scheduled. 

5.2 X-ray data analysis 

Each 4-tile observation with the Swift-XRI is automatically analysed at the UK Swift Science Data 
Centre. The first step of the analysis is to combine the 4 tiles into a single image and to apply the 
source detection algorithm which has been developed for the ISXPS catalogue [36]. It consists of 
filtering the data, creating exposure maps, locating and characterising sources. The source positions 
are generally determined using the on-board star trackers with an accuracy of 3.5 arcsec [37]. How¬ 
ever it is sometimes possible to enhance the astrometry either by using stars in the field of view of 
the UV/optical telescope [38, 39] or by matching the detected XRT sources with the 2MASS cata- 
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Figure 8. Each Swift-XRT observation of an ANTARES trigger consists of 4 tiles (black circles), which covers 
an area of radius of ~ 0.4°. With such a mapping, 72 % of the the bi-directional uncertainty of a TAToO alert 
is covered. 



logue [36]. The position with the smallest error is reported. 

A detection flag is then assigned to each detected source to determine the probability of the source 
being spurious. This flag can be either Good (0.3 % of sources flagged as Good are spurious), Rea¬ 
sonable (7 % of false positive) or Poor (35 % of false positive). More details can be found in [36]. 

Light curves and spectra are also built for each detected source using the tools described by 
[39,40]. A cross-correlation with the X-ray master catalogue, the ISXPS catalogue and SIMBAD^ is 
then performed. The results of this processing, giving details of the newly discovered X-ray source(s) 
and any catalogue matches, are presented in webpages dedicated to each ANTARES alert. 

Finally, two tests are automatically carried out on each uncatalogued X-ray source to determine 
if it could be the counterpart to the ANTARES trigger. The first test relies on brightness and variability 
measurements: if a source is so bright that it should have been discovered previously^, it means that 
the source is assumed to be new and therefore potentially the counterpart. The other possibility is 
if the source shows variability, and in particular significant signs of fading^, it is considered as a 
potential counterpart. The second test is based on the probability that the uncatalogued X-ray source 
is unrelated to the neutrino trigger. The number of expected serendipitous sources with a count-rate 
at least that of the detected source in each 4-tile observation is calculated by combining the minimum 
exposure necessary to detect the XRT source, the sky area, the expected sky density of sources and 
the completeness of the detection system. Details can be found in [42]. 


2 WWW.simbad.u-strasbg.fr 

^The mean XRT count-rate must be more than 1 cr above the Rosat All-Sky Survey (RASS) 3 cr upper limit [41] 
calculated at the location of the X-ray source. The Rosat upper limit is derived using RASS images, exposure maps and 
background maps. 

''The last bin of the source light curve is at least 3 cr below the first bin. 
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5.3 Results 


The Swift-XRT has responded to 7 neutrino triggers from the ANTARES detector between mid 2013 
and the beginning of 2015 (table 4). These 7 alerts have been observed after a mean delay of 3.6 
hours in 4 tiles with exposure from 0.8 to 1.9 ks. With such exposure times, a sensitivity from 6.0 to 
2.6x10“^^ erg cm“^ s“^ can be reached, respectively. 

The tiled analysis found 20 X-ray sources, among which 2 sources were already known (see 
table 5). Although 18 new X-ray sources have been detected, none of them can be clearly associated 
with the neutrino trigger. These uncatalogued sources were not bright enough to exceed the 3 cr 
RASS upper limit. As the XRT failed to observe an X-ray counterpart to a neutrino trigger, this could 
indicate that the observation was made too late or that the source was outside the field of view of 
the XRT, or that the neutrino trigger was not originated by an X-ray source. Upper limits on the flux 
density from a potential X-ray counterpart can be derived. These limits correspond to the sensitivity 
reached for each 4-tile observation (see table 4). 
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Table 2. Details of the 48 neutrinos which triggered the TAToO alert system and for which early optical and/or 
X-ray images were taken. 


Alert name 
(ANTyymmddA/B) 

Trigger time 
(UT) 

Ra 

n 

Dec 

o 

Trigger type'* 

Nhits 

Amp'’ 

Mine*’ 

ANT 100123A 

07:16:45 

73.071 

-2.991 

HE 

34 

339 

4 

ANT 100725A 

01:16:40 

257.224 

-63.188 

HE 

12 

66 

3 

ANT100913A 

21:55:25 

165.142 

41.391 

HE 

46 

277 

5 

ANT100922A 

11:24:23 

43.578 

13.458 

HE 

19 

158 

3 

ANTI 10305A 

12:57:20 

295.517 

-48.193 

HE 

27 

194 

8 

ANTI 10409A 

04:47:50 

129.365 

-40.207 

HE 

21 

253 

8 

ANTI 10531A 

21:45:16 

326.944 

-8.795 

HE 

54 

439 

8 

ANTI 10923A 

08:26:33 

285.194 

10.244 

HE 

22 

236 

7 

ANT110925B 

20:27:38 

70.998 

15.507 

HE 

50 

356 

11 

ANT111008A 

10:45:40 

248.170 

-25.834 

HE 

36 

336 

9 

ANT111019A 

08:56:15 

321.165 

-0.692 

HE 

32 

185 

8 

ANT111019B 

08:56:15 

322.070 

-0.637 

HE 

32 

185 

8 

ANTllllOlA 

12:30:52 

345.917 

-5.698 

HE 

33 

207 

6 

ANTI 11205A 

21:33:11 

146.411 

-32.829 

Dir 

20 

74 

9 

ANT111228A 

09:18:42 

21 ms 

32.323 

Dir 

11 

90 

7 

ANT120102A 

01:24:21 

72.129 

-59.709 

Dir 

11 

54 

8 

ANT120105A 

18:02:53 

228.824 

-26.121 

HE 

26 

184 

7 

ANT120730A 

00:14::33 

178.088 

-40.096 

HE 

135 

489 

11 

ANT120907A 

01:25:10 

220.171 

-11.595 

HE 

77 

245 

10 

ANT120907B 

10:12:41 

344.750 

30.933 

HE 

76 

242 

9 

ANT 120923A 

10:00:16 

318.596 

-51.071 

HE 

71 

128 

4 

ANT121010A 

06:31:01 

52.841 

-29.216 

Dir 

38 

28 

7 

ANT121012A 

04:42:49 

239.512 

-10.226 

Dir 

47 

42 

9 

ANT121027A 

18:47:58 

105.552 

-4.011 

Dir 

39 

51 

7 

ANT121206A 

10:02:42 

62.001 

2.434 

HE 

104 

158 

10 

ANT130210A 

14:00:10 

185.433 

5.897 

Dir 

52 

61 

8 

ANT130722A 

20:38:24 

74.570 

3.411 

VHE 

87 

188 

9 

ANT130724A 

02:07:20 

199.295 

11.938 

HE 

91 

171 

10 

ANT130915A 

06:24:58 

311.314 

-51.603 

VHE 

146 

485 

10 

ANT130927A 

18:22:56 

106.948 

-68.578 

VHE 

81 

258 

10 

ANT130928A 

20:03:18 

121.865 

-3.156 

Dir 

31 

104 

4 

ANT131027A 

12:13:35 

105.124 

6.422 

HE 

74 

190 

7 

ANT131209A 

01:56:24 

197.235 

-13.345 

HE 

87 

163 

9 

ANT131221A 

08:43:36 

138.015 

-23.930 

Dir 

30 

49 

9 

ANT 140123A 

11:53:45 

105.440 

0.811 

VHE* 

80 

401 

8 

ANT140125A 

03:27:08 

53.002 

-49.491 

HE 

80 

224 

8 

ANT 140203A 

06:48:06 

89.731 

-23.394 

Dir 

54 

93 

9 

ANT 140223A 

01:01:20 

43.712 

-10.802 

Dir 

39 

30 

6 

ANT140304A 

05:19:15 

117.930 

-44.944 

HE 

81 

192 

9 

ANT 140309A 

00:40:10 

52.511 

-12.955 

HE 

71 

275 

7 

ANT140311A 

13:16:20 

237.416 

31.826 

VHE 

96 

324 

11 

ANT 140323A 

15:31:01 

150.919 

-27.371 

Dir 

46 

102 

6 

ant140408A 

11:10:09 

177.461 

-3.301 

Dir 

44 

70 

7 

ANT 140505A 

11:54:47 

203.359 

14.450 

HE 

80 

239 

9 

ANT140914A 

05:13:35 

298.424 

2.677 

HE 

109 

210 

9 

ANT141220A 

14:53:59 

71.679 

-45.968 

VHE 

59 

361 

6 

ANT150122A 

07:09:35 

168.551 

-26.119 

Dir 

28 

55 

7 

ANT150129A 

10:20:12 

34.656 

4.863 

VHE 

92 

207 

10 


“ HE: high energy. Dir: directional, VHE: very high energy. With this last trigger, alerts are sent to 
the optical telescope network and to the Swift-XKF. 

Amplitude in photoelectrons. 

Number of ANTARES lines on which the neutrino event has been detected. 

* Early observations by both the optical telescopes and the Swi/t-XRT have only been made for this 
trigger. 
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Table 3. Details of the 42 neutrino alerts for which early optical images have been taken. 


Alert name 
(ANTyymmddA/B) 

Telescope 

Analysed 

images 

Exposure “ 
(sec) 

Delay 

(mag) 

Av** 

(mag) 

Reference‘s 

(yymmdd) 

jpGRB,v f 

reject 

ANT 100123A 

TAROT 

6 

180 

17h47m 

15.3 

0.2 

100211 

0.00 

ANT 100725A 

TAROT 

6 

180 

lml7s 

16.1 

0.3 

100727 

0.50 


ROTSE 

30 

20 

lml5s 

13.1 

0.3 

100731 

0.12 

ANT 100913A 

TAROT 

6 

180 

llh24m 

17.6 

0.0 

100914 

0.06 

ANT 100922A 

ROTSE 

26 

20 

lh08m 

13.6 

0.5 

101202 

0.00 

ANT110305A 

ROTSE 

29 

60 

4hl9m 

15.7 

0.1 

110315 

0.06 

ANT 110409A 

TAROT 

6 

180 

lm08s 

12.6 

5.6 

110423 

0.04 

ANT110531A 

TAROT 

6 

180 

12h34m 

17.6 

0.1 

110628 

0.06 

ANTI 10923A 

TAROT 

7 

180 

9h58m 

12.8 

3.9 

111001 

0.00 

ANT110925B 

TAROT 

6 

180 

2h01m 

15.2 

1.8 

111125 

0.10 


ROTSE 

30 

60 

50m58s 

13.9 

1.8 

110927 

0.00 

ANT 111008A 

TAROT 

5 

180 

12h53m 

14.3 

2.5 

111009 

0.00 

ANT111019A 

ROTSE 

8 

60 

18h22m 

16.7 

0.1 

111029 

0.02 

ANT111019B 

ROTSE 

8 

60 

19h09m 

16.9 

0.1 

111024 

0.02 

ANTI 11101A 

ROTSE 

8 

60 

13h33m 

17.2 

0.1 

111130 

0.02 

ANT 111205A 

TAROT 

6 

180 

10h05m 

18.2 

0.4 

111207 

0.16 

ANT111228A 

TAROT 

6 

180 

7h44m 

17.0 

0.1 

120124 

0.04 


ROTSE 

8 

60 

7h53m 

16.6 

0.1 

111231 

0.04 

ANT 120102A 

TAROT 

4 

180 

lml7s 

17.0 

0.1 

120129 

0.60 

ANT 120105A 

ROTSE 

8 

60 

17h39m 

16.0 

0.4 

120112 

0.02 

ANT120730A 

TAROT 

26 

180 

20s 

16.9 

0.4 

120807 

0.88 

ANT 120907A 

TAROT 

14 

180 

9m53s 

15.9 

0.2 

120909 

0.31 

ANT120907B 

TAROT 

11 

180 

18hl5m 

17.2 

0.2 

121106 

0.02 


ROTSE 

27 

60 

8h28m 

15.9 

0.2 

121004 

0.02 

ANT 120923A 

TAROT 

6 

180 

15h43m 

18.0 

0.1 

121020 

0.03 

ANT121010A 

TAROT 

24 

180 

25s 

18.6 

0.0 

121209 

0.90 

ANT121012A 

TAROT 

6 

180 

19h06m 

16.5 

0.7 

121016 

0.02 

ANT 121027A 

ROTSE 

8 

20 

14h56m 

13.4 

2.6 

121124 

0.00 

ANT121206A 

ROTSE 

27 

60 

27s 

15.6 

1.1 

121210 

0.62 

ANTI 302 lOA 

ROTSE 

8 

60 

14h46m 

16.5 

0.1 

130216 

0.02 

ANT 130724A 

TAROT 

3 

180 

18h04m 

15.9 

0.1 

130729 

0.02 

ANT 130928A 

ROTSE 

8 

60 

13h49m 

15.9 

0.1 

131006 

0.02 

ANT 131027A 

ROTSE 

8 

20 

18hl4m 

15.0 

0.7 

131112 

0.00 

ANT131209A 

TAROT 

6 

180 

lhl4m 

16.3 

0.1 

131215 

0.14 

ANT131221A 

TAROT 

2 

180 

18s 

16.8 

0.5 

131222 

0.83 

ANT140123A 

TAROT 

23 

180 

13h21m 

16.2 

1.3 

140126 

0.02 

ANT140125A 

TAROT 

6 

180 

lhl4m 

18.1 

0.0 

140128 

0.43 

ANT 140203A 

ROTSE 

8 

60 

19h43m 

14.9 

0.1 

140205 

0.00 

ANT 140223A 

TAROT 

3 

180 

17h08m 

15.9 

0.1 

140224 

0.02 


ROTSE 

3 

60 

31m29s 

14.2 

0.1 

140227 

0.02 

ANT140304A 

TAROT 

18 

180 

25s 

17.5 

0.6 

140304 

0.92 

ANT140309A 

TAROT 

16 

180 

24s 

17.1 

0.1 

140309 

0.88 

ANT 140323A 

ROTSE 

8 

60 

14h47m 

16.2 

0.2 

140325 

0.02 

ant 140408A 

TAROT 

6 

180 

16hl Im 

16.7 

0.1 

140410 

0.02 


ROTSE 

8 

60 

19h07m 

16.0 

0.1 

140410 

0.02 

ANT 140505A 

ROTSE 

2 

60 

17hllm 

14.7 

0.1 

140510 

0.00 

ANT140914A 

TAROT 

13 

180 

lm05s 

16.9 

0.5 

140915 

0.62 

ANT 150122A 

TAROT 

8 

180 

17s 

18.6 

0.1 

150122 

0.90 


“ Exposure of each image. 

Delay in hours, minutes and/or seconds between the neutrino trigger and the first image. 

Limiting magnitude of the first image computed at the 5 cr level and corrected for the galactic extinction. 
Galactic extinction from [34]. 

° Date of the reference image used for the subtraction. 

^ Probability to reject the GRB origin hypothesis for this neutrino trigger (see section 6.1). 
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ler of expected serendipitous sources in the field, at least as bright as the detected object. 



6 Limit on a GRB association 


GRBs are the major candidate sources of high-energy neutrinos among the population of fast transient 
sources. In the standard fireball model of GRB [43], internal shocks between faster and slower 
shells of plasma expulsed in a highly relativistic jet would be responsible for the observed prompt 
gamma-ray emission, resulting in photohadronic interactions. Hadronic models of GRBs [44, 45] 
suggest that there could be a significant proton loading in GRB’s jets. If these protons are sufficiently 
accelerated they can interact with the ambient gamma photon field via py mechanisms [46, 47] . The 
photohadronic inferacfions would produce charged pions, subsequenfly decaying info high-energy 
neufrinos. Defails of fhe physical processes faking place in fhe GRB prompf emission are sfill under 
debafe and fhe neufrino expecfafions are highly dependenf on fhe GRB prompf model used. 

According fo fhe sfandard framework of infernal shocks, fhe expected number of defecfable neu¬ 
frinos from fwo luminous bursfs GRB 110918A, GRB130427A and one fypical bursf GRB081008 can 
be esfimafed: see [48] for defails on fhe calculafion mefhod. The presenfed resulfs use fhe predicfed 
neufrino specfra from fhe NeuCosmA model [46, 49] fhaf simulafes fhe complete py mechanism 
(A-resonance, kaon producfion, mulfiple pions decay) in fhe infernal shock GRB model. Sfandard 
values for fhe unknown physical parameters are assumed such as fhe Lorenfz facfor T = 316, fhe 
baryonic loading /p = 10 and fhe fraction of infernal energy given fo fhe magnetic field and elecfrons 
= eg = 0.1. The expected number of prompf neufrinos for each individual GRB is = 3.4 x 10“^, 
ju, = 6.6 X 10-3 and /i, = 2.3 x lO'^ for GRB 110918A, GRB 130427A and GRB081008, respectively. 
So far, no neufrino signal has been observed in coincidence wifh a GRB. Neverfheless, fhe probabilify 
fo observe n neufrinos wifh a given signal can be estimated wifh a Poissonian disfribufion. The 
probabilities fo defecf af leasf one neufrino from fhe mentioned GRBs are P{X > 1 |3.4x lO'^) ~ 3.4 %, 
P{X > 1|6.6 X 10-3) ~ 0.7 % and P{X > 1|2.3 x 10“3) ~ 0.002 %, respectively. Albeif fhese discovery 
probabilities are weak, one should nol forge! fhaf fhey are affected by sfrong fheorefical uncerfainfies. 

The TAToO program offers an opporfunify fo fesf fhe GRB origin for each defected neufrino 
wifh fhe advanfage of nof depending on an underlying GRB model. The probabilify, Pg, to detect a 
serendipitous GRB afterglow lying in the follow-up region for each neutrino alert is given by: 

P,(GRB, y) = Rgrb x p(ag|GRB) x ^ x Tobs (6.1) 

An 

where Rqrb is the average GRB discovery rate (1000/yr) in the entire sky [50] and p(ag|GRB) is 
the probability of detecting the afterglow counterpart of a detected GRB. This probability is 0.4 
and 0.5 for TAROT and ROTSE respectively [51, 52], and 0.95 for Swift [50]. The solid angle, O, 
viewed by the telescopes is given by O = 27r(l - cos{9)) with 0 ~ 1° for TAROT and ROTSE, and 
B ~ 0.2° for the Swift-XRI. Tabs is the exposure time of the observation of each neutrino alert. The 
optical follow-up lasts about 20 minutes, while the X-ray follow-up lasts 2 ks for each tile. Thus, 
Pi(GRB,y) = 1.2 X 10“® and 1.5 x 10“^ for each alert with an optical follow-up with TAROT and 
ROTSE respectively, and Pi(GRB, v) = 1.8 x 10“’ for each alert with X-ray images. Consequently, a 
discovery of a GRB afterglow would be unambiguously associated with the neutrino alert. 

6.1 Optical GRB association 

Because in this study no optical counterpart has been observed in coincidence with the 42 neutrino 
alerts, the probability to reject a GRB associated with each neutrino alert can be directly estimated. 
To do so, a comparison is made between the optical upper limits obtained for each neutrino alert with 
optical afterglow light curves of GRBs detected from 1997 to 2014, as shown in figure 9. 

The cumulative disfribufion funcfions (CDEs) of afferglow magnifudes (corrected by fhe galac- 
fic extinction) are computed af fimes coincidenf wifh fhe firs! follow-up observafion of fhe neufrino 
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Figure 9. Grey lines: Corrected R magnitude as a function of time for 301 GRB afterglows observed from 1997 
to 2014 by optical telescopes. Red & blue dots: upper limits on GRB magnitudes for neutrino alerts observed 
by TAROT and ROTSE respectively. Each point represents the first image of the observation, corresponding 
to an exposure of 180 seconds for TAROT images, and 20 or 60 seconds for ROTSE images. The horizontal 
dashed line corresponds to the maximum sensitivity of the telescopes. 


alerts. Figure 10 shows these CDFs at typical times t = 30 seconds, 5 minutes, 1 hour and 1 day 
after the GRB in the observer frame. The probability, to reject a GRB origin hypothesis for 

each detected neutrino can be directly extracted from the CDFs, accounting for the field of view of 
telescopes^. The listed in table 3 for each alert. 

According to figure 9, sfringenf upper limifs on fhe early aflerglow magnifude a few minufes 
after fhe neufrino frigger can be sef. For fhe fwo neufrino alerfs ANT121010A and ANT150122A, a 
GRB origin is sfrongly rejected since no early aflerglow al a lime 1 < 30 seconds afler fhe bursl has 
ever been observed al weaker magnilude lhan ~ 18. Excepl for Ihe alerl ANT121206A (low limiling 
magnilude), a GRB association is rejected al a belter lhan 80 % confidence level wilhin 1 minute after 
Ihe neufrino frigger. On fhe conlrary, when fhe optical follow-up slarls more lhan a few minutes after 
Ihe frigger, due lo Ihe sensitivity of Ihe telescopes, Ihe derived upper limils do nol conslrainl a GRB 
origin for Ihe detected neulrinos. 

6.2 X-ray GRB association 

As no X-ray counterpart has been identified in coincidence wilh Ihe 7 neufrino alerfs, an analysis 
similar lo fhal in fhe opfical domain is performed wilh fhe 689 X-ray afterglow lighl curves defecled 
by Swift-XRT from 2007 lo 2015 (figure 1 1). A comparison is made belween fhese lighl curves and 
fhe X-ray upper limils oblained for fhe 7 neufrino alerfs. 

^The probabilities extracted from the CDFs are multiplied by the percentage of the ANTARES PSF covered by the 
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Figure 10. Cumulative distribution of afterglow magnitudes for 301 detected GRBs (figure 9). Each line 
corresponds to different times after burst. The vertical dashed line represents the limiting magnitude of the 
optical telescopes. 
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Figure 11. Grey lines; 689 X-ray afterglow fluxes in the energy band from 0.3 to 10 keV detected by the 
Swift-XRT from 2007 to 2015 as a function of time. The upper limits on GRB fluxes for 7 neutrino alerts 
are represented by red triangles. The horizontal dashed line corresponds to the sensitivity reached with a 2 ks 
exposure. 


The CDFs of the X-ray afterglow fluxes are also computed at times coincident with the first 
XRT observation of neutrino alerts. Figure 12 shows these CDFs at t = 1, 2, 4 and 8 hours after 
the GRB in the observer frame. The probabilities to reject a GRB origin of the neutrino events are 


telescopes (see table 1). 
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calculated as previously^ and are listed in table 4. According to figure 11, if each detected neutrino 
comes from a GRB, strong upper limits on the GRB afterglow X-ray flux can be set, Fx ^ 2.5 x 10“^^ 
erg cm“^ s“'. For the alert ANT 130722A, the probability to reject a GRB origin is as high as 71 % 
with an upper limit on the X-ray flux of ~ 2.74 x 10“^^ erg cm“^ s“' only 1.1 hours after the neutrino 
trigger. Clearly, a GRB origin is very unlikely for these neutrino events, since few GRBs are known 
to have weaker X-ray afterglows than these upper limits. 



Figure 12. Cumulative distribution of X-ray afterglow magnitudes for 689 GRBs detected by the Swift-XRT 
since 2007. Each line represents different times after bursts. The vertical dashed line represents the sensitivity 
reached with a 2 ks exposure. 


7 Improvements for TAToO 

The optical and X-ray upper limits derived from the follow-up of 48 neutrino alerts have proved that 
the TAToO system is able to quickly check the GRB origin of the high-energy neutrinos detected by 
ANTARES. Further improvements, such as starting the follow-up as soon as possible after the trigger, 
would enable more definite conclusions on the origin of these high-energy neutrinos. 

7.1 Optical domain 

While small robotic telescopes, such as TAROT and ROTSE, are well-suited to quickly respond to 
ANTARES alerts, their diameters do not permit magnitudes deeper than R ~ 18.5 to be reached. 
This limits the ability of the TAToO program to investigate a GRB origin of the detected neutrinos 
to within only a few minutes after the ANTARES trigger. A one-metre-class robotic telescope is 
crucial to efficiently extend the optical follow-up of ANTARES alerts with stronger constraints on 
the neutrino/GRB association up to a few hours after the neutrino trigger. 

Since 2013, the Zadko telescope has been included in the optical telescope network [53]. The 
Zadko telescope is a one metre fully robotic telescope located at the Gingin observatory in Western 
Australia, covering a field of view of about 0.15 square degrees. This telescope is very interesting 
for the TAToO program as it is located in the area where most of the ANTARES neutrino alerts are 

®For the X-ray follow-up, about 72 % of the ANTARES PSF is covered. 
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immediately visible (figure 5). Thanks to its larger aperture, about 3 times larger than the TAROT and 
ROTSE telescopes, the limiting magnitude of the Zadko telescope with only 60 seconds of exposure 
is 1.4 magnitudes deeper compared to the TAROT telescopes, with 180 seconds of exposure. This 
represents a gain of about a factor of 4 on the optical sensitivity. Due to its limited field of view, 
observations are composed of 7 tiles of 60 seconds exposure each, covering about 85 % of the neutrino 
alert error box. The best response of Zadko to an ANTARES alert achieved so far is about 70 seconds 
after the neutrino trigger. The combination of its reasonable response, its larger aperture and its good 
location will help to better constrain the origin of high-energy neutrinos. 

In addition, five telescopes from the MASTER network [54, 55] have also joined the TAToO 
program since March 2015. They consist of pairs of fully robotic telescopes with a diameter of 0.40 
m covering a field of view of 8 square degrees per pair of telescopes. They are located in Russia 
and South Africa, and significantly increase the coverage of the TAToO telescope network. The first 
alert sendings with MASTER demonstrate the rapid response time to the TAToO alert of about 24 
seconds after the neutrino trigger, which is roughly the same performance as TAROT. Thanks to this 
new partnership with MASTER, the very early optical follow-up (less than a minute) of ANTARES 
alerts will be considerably more efficient at constraining the origin of each neutrino. 

Einally, the addition of an infrared follow-up to TAToO alerts would be very helpful, since this 
part of the electromagnetic spectrum is less subject to dust absorption. In the case of GRBs, this new 
wavelength domain will allow a larger population of GRBs (e.g. extinguished GRB afterglows) to be 
explored, with a correspondingly higher neutrino detection potential. 

7.2 X-ray domain 

In the X-ray domain, the fastest response to the 7 neutrino alerts was about 1 hour. While strong 
constraints on the GRB origin of the high-energy neutrinos detected are set, a faster response or a 
better X-ray sensitivity would almost undoubtedly allow for the rejection (or confirmation) of the 
GRB origin hypothesis. 

With the next generation of neutrino telescopes, such as KM3NeT (see below), the X-ray 
follow-up will be possible with the arrival of new high-energy missions. In the case of the SVOM 
mission (~ 2021) [56], Target-of-Opportunity alerts will be accepted, and X-ray follow-up of neu¬ 
trino alerts by the MXT instrument [57] is under discussion. The sensitivity of the SVOM-MXT 
(3 X 10“'^ erg cm“^ s“' with 1 ks exposure) would allow roughly the same constraints as for the 
Swift-XRT to be reached. Moreover, with the field of view of the MXT (~ 1° x 1°), only one tile will 
be needed, allowing time to be saved for other neutrino alert follow-ups. 

7.3 The future neutrino detector KM3NeT 

KM3NeT^ is the next-generation of neutrino telescope based on the same detection techniques as 
ANTARES in deep seawater. The size of the KM3NeT detector, at several cubic kilometers, will 
increase the sensitivity to neutrino fluxes by about a factor 50 compared to ANTARES. By convolving 
the expected KM3NeT response function with the neutrino flux expected from the NeuCosmA model, 
the number of high-energy neutrinos, fis, expected from GRB110918A and GRB 130427A can be 
estimated. Eor these two GRBs, = 1 and fis - 0.4, respectively, assuming the muon channel only. 
This expectation is clearly encouraging, and the detection of high-energy neutrinos in coincidence 
with GRBs will be achievable with the future KM3NeT detector. Including the faster response of 
the future X-ray telescopes and a better optical/X-ray sensitivity of the TAToO telescope network, a 
factor of ~ 100 in sensitivity would be gained for probing a GRB origin for high-energy neutrinos. 

’ www.km3net.org 
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8 Conclusion 


Optical and X-ray follow-ups of the ANTARES neutrino alerts have been running stably since 2010 
and mid 2013 respectively. Beginning 2015, 150 and 7 alerts have been sent to the robotic telescopes 
TAROT and ROTSE, and to the XRT telescope on board of the Swift satellite. The main advantage of 
the TAToO program is that it is able to send alerts within a few seconds after the neutrino detection 
with a precision better than 0.5° for high-energy neutrinos. Early follow-up has been performed for 42 
alerts with delays as low as 17 seconds between the neutrino detection and the start of the acquisition 
of the first image. The image analysis has not yet discovered any transient sources associated with the 
selected high-energy neutrinos. Upper limits on the magnitude of possible transient sources have been 
derived. Compared to detected GRB afterglow light curves, the very rapid response time has allowed 
stringent constraints on the GRB origin of individual neutrinos to be placed. Even if the response 
time of the XRT follow-up is larger, the early observations of the neutrino alerts have allowed the 
GRB origin for these 7 neutrinos to be excluded with a high confidence level. 
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